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Heterogeneity of solid neutron-star matter: transport 
coefficients and neutrino emissivity 
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ABSTRACT 

Calculations of weak-interaction transition rates and of nuclear formation enthalpies 
show that in isolated neutron stars, the solid phase, above the neutron-drip threshold, 
is amorphous and heterogeneous in nuclear charge. The neutrino emissivities obtained 
are very dependent on the effects of proton shell structure but may be several or- 
ders of magnitude larger than the electron bremsstrahlung neutrino-pair emissivity at 
temperatures ~ 10 9 K. In this phase, electrical and thermal conductivities are much 
smaller than for a homogeneous bcc lattice. In particular, the reduced electrical con- 
ductivity, which is also temperature- independent, must have significant consequences 
for the evolution of high-multipole magnetic fields in neutron stars. 
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1 INTRODUCTION 

The very extensive pulsar and X-ray source observations 
now being made require, for their interpretation, an under- 
standing of the condensed-matter physics of neutron-star 
solid phases. The radiative opacity of a very thin surface 
layer of depth ~ 10 3 cm and matter density p ~ 10 7 g cm~ 3 
largely determines the temperature difference between the 
surface and the interior of the star (see Potekhin, Chabrier 
& Yakovlev 1997, also Potekhin & Yakovlev 2001, for re- 
cent calculations and reviews of earlier work). Atoms in the 
higher- density part of this layer are completely ionized but 
the electron Fermi momentum is less than or of the order of 
1 MeV/c. The whole layer contains no more than 1O~ 1O M0 
and has physical properties which are important only with 
respect to its radiative opacity and, possibly, in connexion 
with the composition of the neutron-star atmosphere. In 
the next layer, with matter densities up to the neutron-drip 
threshold p n d — 4.3 x 10 11 g cm -3 , depth ~ 10 4 cm and mass 
~ 1O~ 5 M0, the electrons form a relativistic Fermi gas whose 
transport coefficients (electrical and thermal conductivities 
a and k) are dependent on the nuclear composition and de- 
gree of order of the solid. This is also true of the neutron-drip 
region with densities above p n d which occupies most of the 
crust volume (depth ~ 10 cm and mass ~ 1O- 2 M ). The 
fraction of the stellar volume concerned here is so large that 
its electrical conductivity must be relevant to the evolution 
of high-multipole components of the magnetic field. Pinning 
of superfluid neutron vortices by nuclei in this region is be- 
lieved to be the origin of pulsar glitch phenomena (Ander- 
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son & Itoh 1975; Ruderman 1976). The neutron-drip region 
contributes almost all the mechanical rigidity of the crust 
and its failure under Maxwell or other stresses is assumed 
to be involved in a number of X-ray emission phenomena, 
for example, the soft gamma repeaters (SGR; Thompson 
& Duncan 1995, 1996) and the persistent emission of the 
anomalous X-ray pulsars (AXP; Thompson et al 2000). 

The canonical picture of the solid, both above and be- 
low p n d, in a neutron star that has not been subject to a long 
period of accretion since formation, is of a homogeneous bcc 
lattice of even-Z nuclei, locally in complete weak-interaction 
equilibrium. Certainly above, and possibly immediately be- 
low, p n d, the nuclei are those with closed proton shells. The 
equation of state below p„d has been derived by extrapola- 
tions of nuclear parameters from experimentally accessible 
regions of neutron excess (see Haensel & Pichon 1994 who 
also summarize earlier work). Pethick & Ravenhall (1995) 
give a general review of solid phase properties and observe 
that, above p n d, there have been two distinct approaches 
to the problem of deriving the equation of state. Micro- 
scopic calculations of the single-particle states for neutrons 
and protons inside a Wigner-Seitz cell were described in the 
classic paper of Negele & Vautherin (1973) and give a com- 
plete description of the system apart from the superfluid en- 
ergy gap. But this approach has not been followed by later 
workers who have adopted a compressible liquid-drop model 
(CLDM) with various Skyrme pseudopotentials (Lattimer et 
al 1985, Douchin & Haensel 2001). Although the CLDM nu- 
clear charge is a continuous variable, the existence of shell 
effects and of proton pairing require that the crust be com- 
posed of successive homogeneous layers of even-Z nuclei. In 
fact, weak-interaction equilibrium cannot be exact owing to 
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the rapid decrease of weak transition rates as the star cools, 
which is caused by the potential barrier present in those 
transitions to oAA-Z nuclei. Flowers & Ruderman (1977) 
noted that, in consequence, at least a small fractional con- 
centration of nuclei with charge deviation AZ = ±2 from the 
homogeneous lattice must be present in metastable equilib- 
rium as point defects. But it has been usual to assume that 
the value of the impurity parameter, 

Q = ^ ai {Z^Z) 2 , (1) 

i 

defined for a distribution of nuclear charges Zi with frac- 
tional concentrations at and mean Z, is, in most con- 
texts, negligibly small (several orders of magnitude less than 
unity). 

The extent of heterogeneity in Z at densities p < p n d 
has been investigated by several authors. Jog & Smith (1982) 
and De Blasio (2000) have examined the structure of the in- 
terface between successive homogeneous layers. The nature 
of the constraints under which equilibrium is defined was 
considered by Jones (1988) who concluded that a distribu- 
tion in Z must be a feature of the state reached as the star 
cools. More recently, calculations of the point defect concen- 
tration have been made by De Blasio & Lazzari (1998). 

Formation enthalpies were calculated by Jones (1999a, 
2001) for a number of point-defects in solid neutron-star 
matter at densities p > p n d- The enthalpies obtained were 
small and it was argued that an amorphous heterogeneous 
solid phase must be formed and should persist as the star 
cools, the melting transition being replaced by a glass tran- 
sition temperature region. 

The conclusion that heterogeneity in Z exists both be- 
low and above p n d does not seem to have been widely ac- 
cepted. Undoubtedly, the assumption of a bec lattice, ho- 
mogeneous in Z, is attractive because it can be simply 
stated and provides a clear basis for calculations, such as 
those of transport coefficients. It is also the case that ar- 
guments about the role of proton shell-structure and the 
approach to weak-interaction equilibrium were made only 
qualitatively by Jones (2001, hereafter Paper I) and were not 
supported by detailed calculation. Nonetheless, even though 
amorphous heterogeneous structures are unattractive owing 
to their greater complexity, if they represent physical reality, 
it is necessary to consider their effect on the stress-response 
of the solid and to define as well as possible the consequent 
degree of uncertainty in calculations of transport coefficients 
and neutrino emissivities. The present paper contains the re- 
sults of those detailed calculations which were absent from 
the previous papers. In Section 2, we consider how the effects 
of proton shell-structure can be included, quantitatively, in 
the compressible liquid-drop model of nuclei by the Strutin- 
ski procedure (see Ring & Schuck 1980) to obtain formation 
enthalpies for nuclei in the interval 20 ^ Z ^ 50. Section 3 
gives estimates of the weak-interaction transition rates be- 
tween these nuclei and, from an initial temperature of 5 x 10 9 
K, which is in the vicinity of the melting temperature or 
glass transition region of the system, describes how the Z- 
distribution evolves with temperature and time. It also gives 
estimates of the very broad range of possible neutrino emis- 
sivities associated with these transitions which have so far 
been neglected in all published calculations of neutron-star 
cooling. 



This paper is concerned with the nature of the crust 
in isolated neutron stars which lack the long period of 
mass- accretion of binary systems. Factors such as accre- 
tion through fall-back at formation are ignored. Therefore, 
it does not consider neutron stars in those binary systems 
where the rate of mass transfer is large enough to replace the 
whole crust, below and above p n d (see Schatz et al 1999). 
The extent to which the results obtained here may be rele- 
vant to such systems is discussed briefly in Section 6. Under 
the physical conditions considered in the present paper, py- 
cnonuclear reactions were examined previously (Jones 2002) 
and were found to have negligible transition rates in the solid 
phase. The reason for this is that the intermediate state 
formed in a solid by fusion of nuclei with charges Z\ and 
Z2 consists of a monovacancy and a point-defect of charge 
Zx + Z 2 . At the highest matter density, 8.8 x 10 13 g cm~ 3 , 
assumed in Sections 2 & 3, for example, this state has a 17 
MeV formation enthalpy. (This assumes that analogues of 
the standard lattice point-defects exist in amorphous solids, 
though they may be short-lived at high temperatures.) 
States with even higher formation enthalpy would result 
from successive fusion reactions, indicating that processes 
such as pinning-induced nuclear rod formation (Mochizuki, 
Oyamatsu & Izuyama 1997) are not significant. However, 
the direct formation of lower-dimensional nuclear structures 
at temperatures of the order of 10 10 K, in a density inter- 
val between the spherical nuclear phase and the continuous 
liquid core of the star, is predicted for many models of nu- 
clear matter (Lorenz, Ravenhall & Pethick 1993; Oyamatsu 
1993; but see also Douchin & Haensel 2000). Lorenz et al 
noted that the geometrical form of these structures would 
allow weak-interaction transitions but gave no estimate of 
the neutrino emissivity. For completeness, a brief calculation 
of this emissivity is given in Section 4. It is less significant 
than that associated with the region of spherical nuclei. 

The density region below p n d is reconsidered in Section 
5 using the binding energy compilation of M0ller, Nix & 
Kratz (1997), but with inconclusive results for values of Q. It 
is possible to state only that values Q « 1 are probable, with 
Q >• 1 in limited regions. This is not too serious a problem 
for most neutron star calculations owing to the limited depth 
of the region below p n d- Values of Q computed for p > p n d, 
the region which occupies most of the crust volume, are given 
in Table 2 and their significance is considered in Section 6. 



2 SHELL EFFECTS IN THE DISTRIBUTION 
OF FORMATION ENTHALPIES 

Formation enthalpies for point-defect structures of impurity 
nuclei with charge Zi in an otherwise homogeneous bec lat- 
tice of charge Z, were obtained in previous work (Paper I) at 
densities above p nd . The method of calculation followed that 
used earlier for monovacancies (Jones 1999b), the essence of 
which was application of the Feynman-Hellmann theorem 
(Slater 1963) to find the lattice displacements in the vicin- 
ity of the defect. Nuclei were described by the compress- 
ible liquid-drop model (CLDM) of Lattimer et al (1985) 
with these authors' Skyrme pseudo-potential for bulk nu- 
clear matter and their expression for the thermodynamic 
potential per unit area of nuclear surface. Lattice-site dis- 
placements in the vicinity of an impurity are determined, 
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principally, by properties of the high-density relativistic elec- 
trons at p > p n d- The Coulomb-electron stress-tensor has 
isotropic components which are between one and two orders 
of magnitude larger than the off-diagonal. Also, the inverse 
of the electron-screening wavenumber is larger than the bcc 
lattice constant. Consequently, lattice-site displacements are 
such that the electron density, averaged over a volume of the 
order of the Wigner-Seitz cell, adjusts to values almost ex- 
actly equal to the mean electron density of the undisturbed 
lattice. We refer to Paper I for more complete discussions 
of these and other features of the formation enthalpy calcu- 
lations. Our primary assumption about the distribution of 
formation enthalpies is based on these considerations. It is 
that the formation enthalpy for a nucleus of charge Zi in 
an amorphous heterogeneous solid of mean charge Z, or in 
a liquid with the same nuclear charge distribution, is sat- 
isfactorily approximated by that calculated for an impurity 
nucleus of charge Zi in a homogeneous lattice of charge Z. 

The formation enthalpies were given in Paper I with 
reference to that for the homogeneous lattice charge. With 
the exclusion of integral multiples of the neutron and elec- 
tron chemical potentials, they were expressed as Hfz = 
C(Z — Z) 2 for charge Z. Values of the constant C are given 
here in Table 1 at several matter densities with, for conve- 
nience, the parameters of the lattices concerned. The values 
of p chosen exclude the region immediately above p n d be- 
cause it represents a relatively small interval of depth in the 
solid crust. The CLDM parameters used by Lattimer et al fit 
the ground states of laboratory even-even nuclei and so in- 
clude the effects of pairing interactions. Thus our expression 
for Hfz neglects shell effects, including the unpaired proton 
in odd-Z nuclei. The investigation by Negele & Vautherin 
(1973) revealed a proton shell structure very similar to that 
of laboratory nuclei. Although shell energy-differences are 
modified in the neutron continuum, the shell ordering (see 
Figure 5 of their paper) is changed only in that ld3/2 and 
lfs/2 respectively, precede 2si/2 and 2p 3 / 2 - We are unaware 
of any published sequence of single-particle energy levels for 
nuclei beyond the neutron-drip threshold apart from those 
contained in that paper. However, Chabanat et al (1998) 
have shown by calculation of two-neutron separation ener- 
gies that shell effects at neutron numbers N = 50, 82 remain 
significant immediately below the neutron-drip threshold. 
There is no doubt that our expression for Hfz ought to be 
modified by shell structure, but the form of these changes 
was considered only qualitatively in Paper I. 

At matter densities p > p n d, there is a neutron con- 
tinuum with chemical potential p n > 0. It is degenerate, 
except for a small volume with p « p nc ;, and is superfluid 
at temperatures < T < T^ n . Nuclei can be most simply 
viewed as bound states of protons embedded in this system, 
charge-neutralized by an almost uniform relativistic electron 
gas. The formation enthalpy differences Hfz+i — Hfz can, 
in principle, be affected by neutron shell-structure because 
the change in Z will be associated with a change in nuclear 
radius rjv and, possibly, a change in the number of neutron 
single-particle states at negative energy. But we shall assume 
that neutron transitions from the continuum just above the 
zero of energy to states just below do not contribute to dis- 
continuities in the formation enthalpy differences with which 
this paper is concerned. Thus we consider the effects of pro- 
ton shell-structure only, and regard the neutrons within the 



nuclear volume as merely a part of the superfluid contin- 
uum, though with increased density and a locally modified 
superfluid energy gap. 

Shell effects in the formation enthalpies of nuclei with 
20 ^ Z ^ 50 are estimated here using the Strutinski pro- 
cedure, following fairly closely the account given by Ring 
& Schuck (1980). If a purely notional single-particle level 
sequence were to be associated with the CLDM approxima- 
tion, its level density would be a monotonic and smoothly 
varying function of Z or of single-particle energy. Real single- 
particle level densities are not like this. Thus the procedure 
starts from a computed level density g(e) and single-particle 
energy sum E sp and generates an averaged level density g(e) 
and energy sum E sp . Hence the formation enthalpy devia- 
tion from its CLDM value at a specific Z is given by 



Hfz = (-Hfz)cldm + E sp — E sp + e v 



(2) 



The average of this expression, for e~ven-Z nuclei, over an 
interval of Z would then be (Hfz)cldm , the underlying 
CLDM formation enthalpy. The correction terms in this ex- 
pression should, of course, be enthalpies because the weak- 
interaction transitions occur at constant pressure. But the 
error in replacing enthalpies by energies here is not large 
and is certainly less serious than that inherent in the ne- 
glect of configuration mixing, to which we shall refer later. 
The first three terms in the right-hand side of equation (2) 
give the shell-corrected formation enthalpy for even-Z nu- 
clei. The remaining term, nonzero only for odd-Z nuclei, is 
the excitation energy of a single proton quasiparticle. The 
shell-model energy is 



eg(e)de, 



(3) 



where A is the proton chemical potential and 

g(e) = J25(e-e i ) (4) 

i 

is the density of single-particle states. The averaged shell- 
model energy is given by the Strutinsky-averaged density of 
states g defined by the averaging procedure 



9(e) = / g(e')f ( — ] d ' ' ( 



It is 



eg(e)de, 



(5) 



(6) 



where the modified chemical potential A remains to be deter- 
mined. It is clearly necessary that successive applications of 
the averaging procedure should leave g unchanged. A suit- 
able class of functions satisfying this condition is formed 
by products of gaussian functions with generalized Laguerre 
polynomials. We assume a specific order of polynomial, 
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and refer to Ring & Schuck for further details and for a 
tabulation of them. An averaged occupation number hi is 
defined for each proton state, 



f(x)dx 



(8) 
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Table 1. Properties of CLDM lattice nuclei in equilibrium with a uniform neutron liquid of number density 
n^. The matter density is p, and Z is here the uniform nuclear charge, a CLDM continuous variable. The 
Fermi wavenumbers are PFn f° r neutrons within the nuclear volume and pp e for the electrons. The radii of 
the nucleus and Wigner-Seitz cell are rjv and r^g, respectively. The lattice Debye and melting temperatures 
are Trj> and T m . The energy gap of the neutron continuum is A*, and C is the formation enthalpy constant. 
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and the constraint 

Y J ^i = Z (9) 

i 

is used to determine the modified chemical potential A. The 
parameter 7 defines the width of the averaging function. It 
is chosen to satisfy, so far as possible, the condition that the 
averaged energy Espd) should be independent of it. 

The scheme requires a set of proton single-particle levels 
e;. As in Paper I, and in order to make use of the values of 
the parameter C calculated in that work, we obtain approx- 
imate sets of homogeneous bec lattice parameters by using 
the CLDM approximation, closely following Lattimer et al 
(1985). These are given in Table 1 and form a basis for our 
calculation of formation enthalpy differences. For each mat- 
ter density, we obtain a set of e, from the work of Negele 
& Vautherin (1973; fig. 5). Values of the width 7 in the in- 
terval 2 ^ 7 ^ 6 MeV have been investigated. The choice 
is to some degree a matter of compromise because the sta- 
tionarity condition is not perfectly satisfied by E sp ('y) at a 
common 7 for all Z in the range 20 ^ Z ^ 50. Thus we have 
adopted a fixed value 7 = 4 MeV for all sets of e; . The first 
three terms in the right-hand side of equation (2) are then 
calculated from the values of the constant C given in Table 
1 and from the procedure described above. The remaining 
term e pq presents some difficulty. It is recognized that its 
value depends on the extent to which the shell concerned is 
filled but for simplicity, we have assumed a constant value 
for all odd-Z nuclei, with e pq = for even-Z nuclei. 

A typical set of Hfz values given by equation (2) is 
shown in Fig. 1, which is for a matter density p = 3.7 x 10 13 
g cm -3 . The enthalpies exclude integral multiples of the 
electron and neutron chemical potentials. There are pro- 
nounced minima in formation enthalpy at the closed-shell 
charges Z = 20, 34, 40, 50. These minima are also present 
at the other densities listed in Table 1 though with rather 
different relative values of formation enthalpy. For the for- 
mation of a thermal equilibrium population of nuclei, or for 
the calculation of weak-interaction transition rates, we re- 
quire only the differences Hfz+\ — Hfz- Also shown are 
(Hfz)cldm values corrected only by the e pq term. Com- 
parison of these two sets confirms that, as suggested in Pa- 
per I, shell effects are a more important source of formation 
enthalpy differences between nuclei of neighbouring Z than 
the CLDM bulk nuclear matter term. 

The fact that the important closed shells are at Z = 
34, 40 is a consequence of our choice, in Paper I, of the Lat- 
timer et al CLDM which produces Z ~ 35 in the density 
interval considered here (see also Pethick & Ravenhall 1995, 
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Figure 1. Formation enthalpies Hpz are given for nuclei 20 ^ 
Z ^ 50 at matter density p = 3.7 X 10 13 g cm -3 . Those repre- 
sented by open circles have been obtained using the Strutinsky 
procedure whereas the solid circles are the (Hpz)cLDM values 
corrected only for the unpaired proton. In both cases, this term 
is e pq = 0.6 MeV. In the former case, proton shell closures at 
Z = 20, 34, 40, 50 are associated with very marked enthalpy min- 
ima. 



Fig. 2). If, instead, the Douchin & Haensel CLDM which 
gives typically Z w 45 had been used, we anticipate that 
the general form of Fig.l would be unchanged except for a 
displacement to higher Z, with the important closed shells 
being those at Z = 40, 50. The shell effects of this paper have 
been obtained using a specific scheme of calculation and it 
is worth considering the extent to which they are physically 
realistic. The values of p concerned are subnuclear so that 
the definition of single-particle levels, as in free (laboratory) 
nuclei, is not unreasonable. The form of single-particle wave 
functions changes completely, owing to spatial quantization, 
in the progression to a new shell and this is reflected in 
the level density discontinuity which produces shell effects. 
But in reality, single-particle or hole states are certainly 
modified by residual nucleon-nucleon interactions which in- 
troduce components of the same angular momentum and 
parity but with more complex particle-hole structure. Our 
expectation is that this configuration mixing would change 
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the relative spacings of the e; so as to reduce the values of 
E sp — E sp generated by this procedure. Consequently, the 
extent of shell structure seen in, for example, Fig. 1, should 
be treated as an upper limit to the true contribution of shell 
structure to formation enthalpy differences between nuclei of 
neighbouring Z. A further reason why these results should 
be regarded as no more than a guide to true formation en- 
thalpy differences is that their Z-average (in Fig. 1) does not 
appear to conform well with (Hfz)cldm, possibly because 
the CLDM parameters are to some extent inconsistent with 
the Negele & Vautherin levels. 



3 WEAK-INTERACTION TRANSITION 
RATES AND COOLING 

Given the formation enthalpies for nuclei in the interval 
20 ^ Z 50, an estimate of the initial condition as the 
star cools can be obtained by assuming that, above a cer- 
tain temperature, weak-interaction transition rates are large 
enough to maintain approximate local thermal equilibrium. 
We assume here, quite arbitrarily, that this temperature is 
To = 5 x 10 9 K. The arbitrary nature of our assumption 
follows from the difficulty in calculating transition rates at 
temperatures T ~ To where neutrino phase-space occupa- 
tion numbers cannot be assumed to be zero. It is also nec- 
essary to make the simplifying approximation of neglecting 
nuclear excited states. Individual nuclear partition functions 
are then equal to 2 J + 1, where the nuclear spin J is de- 
rived entirely from the protons. Unpaired neutrons, or neu- 
tron quasiparticles at T < T cn , are viewed as excitations 
of the neutron continuum rather than of individual nuclei. 
The further evolution of the system at T < To depends on 
the specific heat, the neutrino emissivities, and on the set of 
weak-interaction Z ^ Z + 1 transition rates. 

Equations (A4)-(A9) of Paper I give the Z -> Z + 1 
transition rate from an initial proton closed-shell nucleus. 
A proton is created in a new shell of angular momentum j 
with, in the superfluid case at T < T cn , either the creation 
or annihilation of a neutron quasiparticle. The conservation 
equations for these two energetically distinct processes are 

± e n — Hfz+i — Hfz + e e + £v, (10) 

in which the electron and neutron energies are measured 
from their chemical potentials fi e and ^i n , with e n > A„, 
where A n is the neutron energy gap. These are a form of 
direct Urea transition because the Fourier transform of the 
proton wave function always has a finite amplitude at the 
wavenumber p n — p e p„ necessary for momentum conser- 
vation. The transition rate, summed over all states in the 
new shell, is the product of a rate constant Fo and a phase- 
space integral (equations A7 and A8) . We refer to Paper I for 
further details. In the more general case of a partially filled 
shell, the protons are assumed to be paired into states of zero 
angular momentum. The rate constant Fo is then multiplied 
by a j-dependent factor. For example, in a Z — > Z + 1 tran- 
sition in which the shell initially contains two protons, the 
factor is (1 — 2/ (2j + 1)). A second Z — > Z + 1 case is that in 
which the shell initially contains an odd number of protons. 
For the case of a single proton, the factor is 2/(2j + l) 2 . 

Nuclear spins obtained directly from the Negele & Vau- 
therin shell ordering are assumed in the present paper, but 



with some reservation. The sequences of ground-state spins 
and parities tabulated by M0ller et al (1997) for nuclei close 
to neutron instability, and hence relevant below neutron- 
drip, are not at all simple, possibly owing to the significant 
nuclear deformation. It is not obvious that, above neutron- 
drip, a system of Z protons embedded in a neutron contin- 
uum will be without deformation and consequent complex- 
ities. But the conclusions which will be described in this 
Section concerning Z-heterogeneity are so clear that they 
would not be affected by changes in the detail of shell or- 
dering. 

The parameters used in the evaluation of transition 
rates at a given temperature are those listed in Table 1 ex- 
cept that we assume different neutron energy gaps and crit- 
ical temperatures in the nuclear interior and in the contin- 
uum. The neutron effective mass is m* = 0.8m n . Following 
a local density approximation (LDA), values of the neutron 
Fermi wavenumber pFn and the zero-temperature neutron 
energy gap A n are those for an infinite system at the den- 
sity of the nuclear interior. The critical temperature T cn is 
that for an isotropic BCS superfluid. The finite-temperature 
energy gap needed for calculation of the phase-space inte- 
gral is conveniently obtained from tabulated results given by 
Rickayzcn (1965), as is the BCS neutron specific heat. The 
specific heat of the system has components for the electron 
gas, the solid, and for the neutrons (normal or BCS super- 
fluid). Although the solid is amorphous and heterogeneous, 
it is represented by a simple Debye model with a single tem- 
perature To- 

The neutron energy gap in neutron-star matter is not 
well known. Different methods of calculation lead to a wide 
range of values, (see the reviews of Pethick & Ravenhall 
1995, Heiselberg & Hjorth- Jensen 2000). However, methods 
that introduce corrections such as medium polarization to 
the calculation inevitably seem to suppress the energy gap 
as shown, for example, by Shen et al (2003) who also give a 
brief review of recent work on this topic. For this reason, we 
have assumed that the energy gap obtained by Ainsworth, 
Wambach & Pines (1989) applies to the neutron continuum 
A^. This affects only the specific heat. The neutron energy 
gap effective in the nuclear interior and the proton quasipar- 
ticle energy are even more difficult to estimate in neutron 
star matter. The neutron coherence length is at least of the 
same order of magnitude as the nuclear radius, so that the 
proximity effect is certain to be significant for the neutron 
gap, as has been emphasized by Pethick & Ravenhall. Hence 
the value of the LDA assumption must be limited. These 
authors also observe that pairing in the interior of ordinary 
nuclei is a kind of proximity effect in that attractive pairing 
in the relatively large, low-density, surface volume influences 
pairing in the whole. The absence of this low nucleon- density 
region owing to the presence of the neutron continuum may 
then affect the proton pairing e pq . The uncertainties of pre- 
diction are such we think it necessary to treat the interior 
neutron energy gap and the proton quasiparticle energy as 
unknown parameters which are both assumed to be within 
the interval < e pq , A n < 1 MeV. 

Routines for the evaluation, at any temperature, of the 
specific heat and of all transition rates for 20 ^ Z ^ 50 
with the proton level sequence lf 7 /2, H5/2, 2p 3 / 2 , 2pi /2 , lg 9 / 2 
allow the change in ^-distribution to be followed as the 
star cools. They also give the neutrino emissivity. We have 
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Table 2. Values of the impurity parameter Q and mean charge 
Z: Q = Q sc with shell corrections; Q = Qnsc with no shell 
corrections. 



6 8 10 

log t (s) 

Figure 2. The left-hand set of curves shows the cooling as a 
function of time at a matter density of 3.7 X 10 13 g cm -3 for for- 
mation enthalpies given by the solid circles of Fig. 1. The cooling 
is adiabatic except for the neutrino and antineutrino emissivities 
described by equation (10). Neutrino pair production by electron 
bremsstrahlung and by neutron quasiparticle annihilation are ne- 
glected here. Allowance for these processes would produce faster 
cooling, particularly at t ~ 10 10 s. The right-hand curves measure 
the movement of nuclear charge (the flux f w i per unit nucleus) 
from Z = 40 toward the closed shell at Z = 34 as a consequence 
of the weak interaction. It is a function of temperature, but is 
represented in the figure as a function of time by means of the 
left-hand set of curves. The product tf w i is always some orders 
of magnitude smaller than unity, showing that the movement of 
charge during cooling is insignificant. The solid, large and small- 
dashed curves are, respectively, for the parameter sets: A n = 0.2, 
e vq = 0.4 MeV; A n = 0.4, e vq = 0.4 MeV; A n = 0.4, n vq = 0.6 
MeV. 



adopted an arbitrary fixed initial temperature To — 5 x 10 9 
K, which is close to the bcc lattice melting temperatures 
given in Table 1. The cooling curves for a matter density 
of 3.7 x 10 13 g cm~ 3 (adiabatic except for neutrino and 
antineutrino emission through the processes described by 
equation 10) are shown in Fig. 2. for the set of Hfz val- 
ues given in Fig.l which have no shell correction except for 
the e pq unpaired-proton term. Three sets of the parame- 
ters t pq and A n are considered. The early parts of these 
curves are not entirely reliable because our calculations ne- 
glect neutrino opacity and so overestimate transition rates 
and emissivities at temperatures significantly above 10 9 K. 
The right-hand scale gives a measure of how rapidly, at any 
given temperature, the weak interaction is changing the nu- 
clear Z-distribution by moving nuclei between two of the 
minima in Z which appear in Fig. 1. It shows the flux f w i 
between Z — 40 and Z = 34, calculated at Z = 37 and 
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normalized per unit nucleus in the system, for the cooling 
conditions assumed here. This quantity is a function of tem- 
perature, apart from a relatively small dependence on the 
changing Z-distribution. It is given here as a function of 
time by using the left-hand sets of curves in the Figure. Us- 
ing the (right-hand) sets of curves so obtained, it can be 
seen that the product tfwi is always some orders of mag- 
nitude smaller than unity. (Fig. 2 allows recovery of f w i as 
a function of T so that it is possible to estimate the time 
required for weak-interaction equilibrium at a fixed temper- 
ature.) Weak- interaction rates are so low that even with the 
extreme assumption of no shell corrections, the distributions 
of fractional concentrations for even-Z nuclei change very lit- 
tle with time provided e pq ^ 0.4 MeV. In the limit of large 
t, the resulting values of the impurity parameter Q are not 
much smaller than those (Q = Q m in Table 2 of Paper I) 
for thermal equilibrium at the melting temperatures without 
shell corrections and with e pq = 0. They are shown here in 
Table 2. In the shell correction cases, there is naturally some 
movement of Z toward the magic numbers but the resulting 
values of Q are independent of A n and e pq to such an extent 
that, given the very approximate nature of our calculations, 
it is not useful to include in the tabulation the values of these 
two parameters used, which are those of Figs. 4-6. The same 
comment can be made for the no shell-correction cases pro- 
vided, as we stated previously, that t vq ^ 0.4 MeV. Lower 
values of this parameter, in the special case of the absence of 
shell corrections, produce potential barriers small enough to 
allow movement of Z toward Z . But we do not attempt to 
give Q-values for this set of assumptions which are very un- 
likely to be physically realized. The Q-values of Table 2 and 
their significance will be considered further in Section 6. But 
it is not realistic to neglect shell effects completely. These 
must be present although, owing to uncertainties such as 
configurational mixing (Section 2), the formation enthalpy 
differences they introduce may be smaller than those shown 
in Fig.l. Thus the potential barriers in the weak-interaction 
paths between the closed shells at Z = 20, 28, 34, 40, 50 must 
be much higher than assumed in our calculations of f W i , and 
the actual values of f m i many orders of magnitude smaller 
as a consequence. 

Sets of curves similar to those of Fig. 2 exist for the 
other densities listed in Table 1. But they are so similar 
that they are not shown here. Our conclusion is that, even 
with allowance for the uncertainties inherent in the shell- 
effect formation enthalpy differences calculated here, solid 
neutron-star matter cools to a system heterogeneous in Z. 
Other sources of error are much less significant. For ex- 
ample, the Bessel function bound-proton states assumed 
in equations (A4)-(A6) of Paper I must overestimate the 
true transition rates, as would our neglect of neutrino opac- 
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Figure 3. The specific heat per average Wigner-Seitz cell is 
shown as a function of temperature, in units of kg. The curves 
labelled 1 — 3 are, respectively, for matter densities 1.6,3.7, and 
8.8 X 10 13 g cm -3 . The neutron continuum energy gaps A^ are 
those given in Table 1. The internal gap A n = 0.2 MeV. 



ity at T > 10 K. But modest errors are not important 
because these transition rates determine both the cooling 
and change in Z-distribution of the system. Very broadly, 
we can see that a sufficient condition for maintaining Z- 
heterogeneity during cooling is that the thermal energy per 
average Wigner-Seitz cell should be smaller than the mean 
energy, of the order of 6ksT, removed by neutrinos or an- 
tineutrinos in a weak transition. The specific heat per av- 
erage Wigner-Seitz cell is shown in Fig. 3 as a function of 
temperature. It demonstrates that the thermal energy sat- 
isfies this condition easily at T ~ 10 9 K but not at T> 10 9 
K (the neutron contribution is large at T ~ T c e n ). It is also 
true that we have neglected other well-established sources of 
neutrino emissivity which, if included, would cool the system 
more quickly and so assist in maintaining heterogeneity. 

The neutrino emissivities are shown in Figs. 4-6 for the 
matter densities listed in Table 1. In each case, they are 
given both with and without the Strutinsky shell corrections 
obtained here and for sets of values of t pq and A„. The prin- 
cipal well-established neutrino-emissivity processes, for the 
inner crust, in studies of neutron-star cooling are neutrino- 
pair production by electron bremsstrahlung and by neutron 
quasiparticle annihilation in the superfluid. We emphasize 
that these emissivities have not been included in our cool- 
ing calculations but, for reference purposes, they are shown 
separately, obtained from Fig. 4 of the paper by Kaminker, 
Yakovlev & Gnedin (2002). Obviously, the potential barriers 
which separate formation enthalpy minima at closed shells 
are of crucial importance in determining emissivities. It is 




8.4 8.6 
log T (K) 

Figure 4. Neutrino emissivities at a matter density of 1.6 X 10 13 
g cm -3 are shown as functions of temperature for the weak- 
interaction processes described by equation (10). The formation 
enthalpies are obtained by the inclusion of shell corrections (lower 
set of curves) and without these corrections (upper set of curves). 
Neutrino-pair production by electron bremsstrahlung and neu- 
tron quasiparticle annihilation are excluded from these curves, 
but their emissivity is shown separately for reference purposes. 
The curves labelled a — c are, repectively, for the following pa- 
rameter sets: An = 0.2, e pq = 0.4 MeV; A n = 0.4, e pq = 0.4 
MeV; A n = 0.4, e pq = 0.6 MeV. 



possible to say with confidence only that the emissivities 
shown are likely to form upper and lower bounds for the true 
values. As they differ by at least several orders of magnitude, 
this is not a very strong or practically useful statement. But 
given that the solid inner crust is heterogeneous in Z, the 
existence of this source of uncertainty must be accepted. Its 
effect on the surface temperature of the star is difficult to 
estimate without complete cooling calculations. 



4 LOW-DIMENSIONAL STRUCTURES 

Phases in which protons are confined to one or two- 
dimensional structures have been studied extensively and 
their existence in a substantial density interval between the 
spherical nuclear phase and the liquid core depends on the 
form of the Skyrme pseudopotential assumed in deriving the 
equation of state. Both Lorenz et al (1993), using the inter- 
action of Lattimer et al, and Oyamatsu (1993) find such 
structures within a significant density interval. A survey of 
work on this topic has been given by Pethick & Ravenhall 
(1995). But more recent calculations by Douchin & Haensel 
(2000), using a different Skyrme interaction, find only the 
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Figure 5. Neutrino emissivities at a matter density of 3.7 X 10 
g cm -3 are shown with labelling and parameter sets as for Fig 
4. 



phase of spherical nuclei. Investigations since then (Watan- 
abe & Iida 2003) have shown that the inclusion of electron 
screening tends to increase the density interval occupied by 
any such phase. The paper by Lorenz et al notes that these 
structures give rise to new weak-interaction processes but 
gives no further details. These are no more than a differ- 
ent form of the direct Urea transition considered in Section 
3. For completeness, the neutrino emissivity has been cal- 
culated here for one of these cases, allowing a comparison 
with the range of emissivities obtained for the spherical nu- 
clear phase. The one-dimensional system considered is that 
in which the Wigner-Seitz cell is an infinite slab of thickness 
2dws containing neutrons and relativistic electrons. Protons 
are confined within a slab of thickness 2djv symmetrically 
positioned inside this cell. Thus the component of the pro- 
ton wave function in the variable perpendicular to the slab 
has Fourier transform ip m with quantum number m and fi- 
nite amplitude at the wavenumber necessary for momentum 
conservation. We assume that, with this limitation, neutrons 
and protons each form an isotropic BCS superfluid, with en- 
ergy gaps A„ |P . The neutrino emissivity per unit volume is 
then 

(1 + 3Ci)G 2 cos 2 9 C f , 3 , 2 l3 , 3 

A > 1 d i p n d 2 qd i p e d A p„ 



(2-ir) 11 MwS 

^ \^m(Pn± -Pe±)f 5 {2) (p n - p e - q) 

m 

y ) j e„?4«p7i, ra (l - ra p )(l - n e ) 
S(e„ — e p — e e — e v ). 
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Figure 6. Neutrino emissivities at a matter density of 8.8 X 10 
g cm -3 are shown with labelling and parameter sets as for Fig. 
4. The additional set of broken curves are the neutrino emissiv- 
ity L BV given by equation (11) for a one-dimensional slab nuclear 
structure. The curves labelled d — f are, respectively, for the pa- 



rameter sets: A n = 
MeV; An = 0.4, A K 



0.2, A p 
= 0.6 MeV. 



0.4 MeV; A n = 0.4, A„ 



0.4 



(11) 



In this expression, G is the muon decay constant, 9 C is the 
Cabibbo angle, and Ca = 1.25 is the ratio of axial vector 
to vector coupling constants. The wavevector q lies in the 
plane of the slab. The neutrino or antineutrino energy is 
t v . Quasiparticle or electron occupation numbers are n n ,p,e 
and £n,p,e are the energies, referred to the chemical poten- 
tial. The Bogoliubov coefficients for the proton and neutron 
quasiparticle states are u n and v p . The terms shown in this 
expression are for the process of neutron quasiparticle an- 
nihilation with proton quasiparticle and electron creation. 
However, the unlabelled summation states that the set of 
8 terms, derivable from the creation or annihilation of the 
electron and of the neutron and proton quasiparticles, are 
included. 

Numerical evaluation of this somewhat untidy expres- 
sion gives the emissivities shown in Fig. 6 for a matter den- 
sity of 1.59 x 10 14 g cm" 3 (djv = 3.62 fm, d ws = 10.1 fm) 
for several sets of values of A n ,p. Values of the emissivity do 
not vary greatly with matter density in this phase or with 
the parameter djv which is itself very much dependent on 
CLDM details. Emissivity calculations with different values 
of djv show that uncertainties arising here are small com- 
pared with those caused by our lack of knowledge of the en- 
ergy gaps A n , p . In general, the emissivities are rather smaller 
than those found in the spherical nuclear phase, the reason 



Heterogeneity of solid neutron-star matter 9 



being that the Fourier transforms of the proton states in 
that phase have broader distributions of momentum. 



5 THE OUTER CRUST 

The canonical crust structure at densities p < p n d is ob- 
tained by minimization of the Gibbs free energy per nu- 
cleon. It consists of successive spherical shells, each of a bcc 
lattice homogeneous in the nuclear charge Z neutralized by 
a relativistic electron gas. The most recent calculations are 
those of Haensel & Pichon (1994) who also review earlier 
work. The zero-temperature structure of the interfaces be- 
tween these shells was examined by Jog & Smith (1982) and 
found to consist of extremely thin layers of interpenetrating 
simple cubic lattices of the two nuclear charges concerned. 
Formation, however, must be considered at the lattice melt- 
ing temperature T m , or above. The reduction in free energy, 
derived from the melting temperature entropy of mixing for 
the two nuclear species, then gives interface layers of signif- 
icant thickness in which one nuclear species is present as an 
impurity in a lattice formed by the other (De Blasio 2000). 

Consideration of the later stage of nuclear formation at 
finite temperature gives quite separate grounds for excluding 
homogeneous lattices (Jones 1988). At temperatures in the 
vicinity of T m , nuclei are in state of partial thermal equilib- 
rium with the electrons and with a low-density Boltzmann 
gas of neutrons. (Weak-interaction transition rates, even at 
T m , are too small to guarantee equilibrium Z- values, but the 
transition rates for the emission or absorption of neutrons 
are high and so maintain strong-interaction equilibrium with 
the Boltzmann gas at any density.) From the definition of 
the neutron-drip density p n£ j, the neutron chemical potential 
in this region is p n ~ Pn < 0, referred to the rest energy as 
zero. The equilibrium number density of the Boltzmann gas 
is 

2 

An 



n n = tj exp(^), 



(12) 



where 
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(13) 



and = UbT. Thus it decreases rapidly as the star cools. 
But the nature of the system imposes an additional con- 
straint on the Gibbs function minimization. The neutron 
scattering cross-section is so large that, once nuclei have 
formed, diffusion over macroscopic distances of the order of 
the shell depth is not possible within the short time permit- 
ted by neutron condensation on to nuclei. This occurs at a 
temperature T > To, where To is the temperature (not well- 
defined) at which weak-interaction equilibrium fails. Conse- 
quently, it is necessary to impose the further condition, at 
temperatures T ~ T m , that the mean nucleon number Aws, 
per Wigner-Seitz cell, should be constant. We refer to Jones 
(1988) for further details. 

The method of calculating equilibrium nuclear number 
densities differs from that of Paper I which was for densi- 
ties above p„d- The pressure at p < p„d is small in nuclear 
structure terms, being almost entirely that of the degener- 
ate electron gas. Nuclei in this region have definite A, Z and 
can be assumed to have binding energies B(A,Z) identical 



with those of the free (terrestrial) state. The minimization 
of the Gibbs function G for a fixed number of baryons is at 
constant p e and hence approximately at constant pressure. 
It is, 
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(14) 



where mj is the number of nuclei with mass number Ai and 
charge Zj, and A is a Lagrange multiplier. The equilibrium 
condition for the nuclear chemical potential fuj is, 



Aij + {Aws - Ai)(p" n + k B T) + Z 3 p e + A = 0. 



(15) 



(Without the constant Aws constraint, the nuclear chem- 
ical potential would be p,ij — Aipf n — Zjp e .) The neutron 
chemical potential is not an independent variable in equa- 
tion (15): it is a function of Aws through the relation 
Aws — A = n%Vws, in which A is the mean nuclear mass 
number and Vws the Wigner-Seitz cell volume. To a satis- 
factory approximation, the nuclei are independent systems 
and their equilibrium number densities are given by an ex- 
pression entirely analogous with equation (12), 

2, 

riA,z = Y3 — ex P (Pifiii + z ( m n - m p ) + B - E W s)) , (16) 

A A,Z 

where Z is the nuclear partition function (normalized so that 
at T = it is 2 J + 1), B is the nuclear ground-state binding 
energy, and Ews is the Coulomb energy of the Wigner-Seitz 
unit cell. The presence of A allows two finite number den- 
sities in the zero-temperature limit, as is necessary because 
Aws is in general not an integer. Evaluation of equation (16) 
at a standard temperature To , close to T m , for even-even 
nuclei with Z — 1 and binding energies from the recent com- 
pilation of M0ller et al (1997) gives estimates of the mass 
number and charge heterogeneity expected at densities be- 
low p nd . 

It is unfortunate that the results are inconclusive. The 
source of the problem is that the values of Aws are deter- 
mined at those temperatures T > To existing at the forma- 
tion of nuclei, which are high and poorly known. The bulk 
transport of neutrons over macroscopic distances, of the or- 
der of shell depths, is not possible in the short time which 
elapses before rapid cooling produces almost complete con- 
densation through nuclear capture (Jones 1988). It has to 
be accepted that the values of Aws, constant at T < T, are 
virtually unknown. Evaluations of equation (16) for matter 
densities p ~ 1 — 3 x 10 11 g cm" 3 , temperatures To = 5 x 10 9 
and 10 10 K, and a wide range of values of Aws — A give a 
common picture. As Aws ~ A increases, the equilibrium nu- 
clei change from a group associated with the neutron A^ = 50 
closed shell to a group near N = 82. The neutron chemical 
potentials at changeover vary from —1.3 MeV at 3 x 10 11 g 
cm" 3 to -3.6 MeV at 1 x 10 11 g cm" 3 , for T = 5 x 10 9 
K. The corresponding values at 10 10 K are —1.8 and —4.1 
Mev, respectively. The changeover from N = 50 to N = 82 
is similar to that found by Haensel & Pichon (1994) in their 
study of zero-temperature equilibrium. The problem is that 
the Aws constraint which was not considered by them does 
not permit us to estimate where it occurs. This complexity 
is unfortunate, but we believe that it is real. Qualitatively, 
the consequences for Q-values are as follows. The N — 50 
region has Q < 1, and the A^ = 82 region, Q ~ 1. Values 
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Q> 1 exist in the changeover region owing to the very large 
differences in Z which are present there, as found earlier by 
De Blasio (2000). Detailed calculations of weak-interaction 
transition rates confirming the absence, at T < To, of weak- 
interaction equilibrium have not been made for p < p n d, 
but the qualitative criterion given in Section 3 is satisfied. 
In particular, values Q>1, where they exist, will certainly 
remain frozen in. 



6 TRANSPORT COEFFICIENTS AND 
CONCLUSIONS 

The main conclusion of this paper is that nuclear-charge het- 
erogeneity exists in the solid phase of isolated neutron stars. 
There must be many reservations about the procedures de- 
scribed in Sections 2 & 3 which have been used to obtain 
this result. We attempt to summarize them here and in each 
case give reasons why they should not be viewed as seri- 
ous. Firstly, the shell-corrected formation enthalpy values, 
Hfz, are based on the Negele & Vautherin shell ordering 
and spacings. Nuclei close to neutron-instability are known 
to be deformed, with ground states of some complexity (see 
M0ller et al 1997). It is quite possible that these features 
are also present above pnd- Therefore, the values shown in 
Fig. 1 may not give a true picture of the potential barriers 
which slow weak-interaction transitions. But it is the case 
that even quite small barriers slow weak transitions suffi- 
ciently. The transition rate calculations, following the pro- 
cedure of Paper I, are based on elementary single-particle 
proton wave functions and neglect Coulomb corrections to 
the electron function, which are only moderate at the val- 
ues of fj, e considered. But here, the effect of more complex 
ground states would inevitably be reductions in transition 
rates which would be largely neutral in effect because the 
rate of cooling by neutrino emission is also reduced. The re- 
sults, of which those in Fig. 2 are an example, show that 
the failure of weak-interaction equilibrium during cooling at 
T < To is so clear that the above uncertainties are not sig- 
nificant. The choice of To = 5 x 10 9 K as the lowest temper- 
ature at which complete equilibrium remains is, of course, 
arbitrary. There is no doubt that a higher value would give 
more heterogeneity in Z, with larger values of Q, and that 
calculations made by the present methods would show a fail- 
ure of weak-interaction equilibrium at temperatures below 
it. But at these very high temperatures, the neutrino opac- 
ity would cease to be negligible, as assumed in Section 3, 
so that the cooling rates shown in Fig. 2 would be over- 
estimates. Obviously, the Q-values given in Table 2, under 
different stated assumptions, for each of the matter densities 
listed in Table 1, should be seen as little better than order 
of magnitude estimates. But these satisfy Q > 1 in a way 
which does not depend systematically on the details of shell 
effects assumed. They are some orders of magnitude larger 
than the values Q <C 1 previously assumed in the standard 
view of solid neutron-star matter. Although the relevant cal- 
culations have not been made, we suggest there is no reason 
to think that this conclusion would be changed by replacing 
the Lattimer et al CLDM with that of Douchin & Haensel 
(2001). In the latter case, the important proton closed shells 
would be those at Z — 40, 50 but we anticipate that the 



orders of magnitude of the Q-parameter obtained would be 
the same. 

The procedures of Sections 2 & 3 are concerned merely 
with establishing Z-heterogeneity. The small-scale structure 
of the solid is a different question. It is not obvious that dif- 
fusion at temperatures in the glass transition region, which 
we assume to be near the homogeneous lattice melting tem- 
perature Tm, could produce localized order extending over 
linear dimensions much greater than lO 1 ^ 2 inter-nuclear sep- 
arations. The results of rudimentary formation enthalpy cal- 
culations given in Paper I do not indicate that large chemical 
potential gradients exist to drive such diffusion. Moreover, 
the classical entropy of disorder remains large. For these rea- 
sons, it is believed that although the structure may not be 
exactly that of an amorphous heterogeneous solid, any local 
order will be at most of small linear dimension, analogous 
with nanostructures in ordinary amorphous matter. 

These conclusions can be compared with those of re- 
cent work on neutron stars in binary systems with high mass 
transfer rates, ~ 1O~ S M yr _1 (Schatz et al 1999, see also 
Brown & Bildsten 1998). Stable burning of hydrogen and 
helium occurs near the surface, the most important process 
being rapid proton capture, and is almost complete at a 
depth of the order of 10 s g cm -3 (see Schatz et al, Fig.l). 
The end products have charge distributions which depend on 
the mass-transfer rate, but are typically broad, with Z < 40. 
Continued mass-accretion forces these nuclei to higher mat- 
ter densities, eventually to p > p n d- Schatz et al contrast this 
Z-heterogeneous matter with the solid phase of an isolated 
(primordial) neutron star having the properties described by 
Pethick & Ravenhall (1995). Our conclusion is that the solid 
phases of neutron stars in these different environments are 
broadly similar, with impurity parameters Q ~ 10. It is very 
unlikely that the temperatures of the burning processes will 
be high enough to produce weak-interaction equilibrium at 
matter densities above pnd- The flux f W i, representing the 
movement of nuclear charge toward a closed proton shell, 
is shown in Fig. 2 for the case of no shell correction, apart 
from the proton pairing term e pq . This extreme case is not 
realistic because there is almost certainly some shell correc- 
tion, though perhaps not so large as that shown in Fig. 1. 
The increased potential barriers introduced by the correc- 
tion would reduce values of f wi by many orders of magni- 
tude, to an extent that Z-heterogeneity would be largely 
unaffected. 

It is obvious that the amorphous and heterogeneous na- 
ture of the solid phase of neutron-star matter has significant 
consequences for its mechanical properties and, above p n d, 
for its interaction with supcrfluid neutron vortices. Mechan- 
ical properties will not be considered here except to note 
that the response to stress is not strictly that of an ordi- 
nary amorphous solid because neutron-star matter is very 
far from being absolutely stable and therefore cannot ex- 
hibit brittle fracture (Jones 2003). Superfluid neutron vor- 
tices interact with nuclei through elementary pinning forces 
whose calculation is associated with unresolved problems 
(for a brief partial review, see Jones 2002). It is currently 
not possible to say that the signs or orders of magnitude 
are known with any degree of certainty at a given mat- 
ter density. Both vortex pinning and the dissipative force 
acting on a moving vortex are strongly influenced by the 
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structure of the solid phase. Apart from the kind of low- 
dimensional structure considered in Section 4, very small 
dissipative forces would be possible only for motion through 
a system of large single crystals with very low concentra- 
tions of dislocations and point-defects. The confirmation of 
the amorphous and heterogeneous nature of the solid phase 
rules out this possibility. 

The neutrino emissivity from nuclear weak interactions 
at matter densities below p n£ j has not been calculated ow- 
ing to the difficulty in knowing the correct transition rates 
between nuclei which are close to being neutron-unstable. 
The assumption of superallowed Fermi transitions gives an 
order of magnitude of 10 22 erg cm~ 3 s _1 at 10 9 K for matter 
densities 1 — 3 x 10 11 g cm~ 3 , but allowance for the forbid- 
den nature of the transitions would probably reduce this to 
an emissivity not much different from that for neutrino-pair 
production by electron bremsstrahlung under these condi- 
tions (10 18 erg cm~ 3 s _1 ; see, for example, Gnedin et al 
2001). Neutrino emissivity above p n d is shown in Figs. 4-6. 
Weak-interaction transitions in this region involve the cre- 
ation or annihilation of quasiparticles in the neutron con- 
tinuum. Rate calculations are more straightforward in this 
case, but owing to the uncertainties in formation enthalpy 
differences considered above and in Section 2, the emissivi- 
ties shown in the two sets of curves probably form upper and 
lower bounds for the true values. However, at temperatures 
T ~ 4 x 10 8 K, they can be larger than those previously 
assumed. The degree of uncertainty here is unfortunate but 
it does seem right that its existence should be recognized in 
neutron star cooling calculations. Neutrino emissivities for 
one-dimensional layer structures (Section 4) are also sub- 
ject to some uncertainty, but are rather smaller than for the 
spherical nucleus phase considered in Section 3. Emissivities 
for a two-dimensional phase have not been calculated, but 
are probably of a similar order of magnitude to those for one 
dimension. 

Electron scattering in a ^-heterogeneous system con- 
tributes electrical and thermal resistivities which are related 
by the Wiedemann-Franz law and are proportional to Q. 
The corresponding electrical conductivity (for zero magnetic 
flux density) is, for relativistic electrons at p ~ 10 7 g cm -3 , 
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(17) 



2.03 in the ultra-relativistic 



in which the parameter Aj n 
case (see Urpin & Yakovlev 1980, Itoh & Kohyama 1993) 
and Z is the mean nuclear charge. For the values of Q ob- 
tained here, the relaxation time underlying equation (17) is 
so short that its product with the electron cyclotron angular 
frequency exceeds unity only for magnetic fields B ~ 10 
G. Therefore, there is no need to distinguish between longi- 
tudinal and transverse components of the conductivity ten- 
sor, except at B > 10 13 G. The conductivity is also classi- 
cal rather than quantum in nature at densities of the order 
of p n d or above (see Potekhin 1999). The electron thermal 
conductivity m is given, in terms of by the Wiedemann- 
Franz law. These transport coefficients are of importance, in 
most cases, at very different times during neutron-star cool- 
ing. The thermal conductivity k is usually of interest at early 
times before the interior of the star, defined for this purpose 
as matter densities p ~ 10 g cm , becomes approximately 



isothermal. The Q-dependent resistivity obtained here de- 
creases k, but the significance of the change is not immedi- 
ately obvious and would require complete neutron-star cool- 
ing calculations for its investigation. At later times, internal 
temperature gradients are very small and it is possible that 
the consequent changes in them would be relatively unim- 
portant. 

The electrical conductivity a is obtained by combin- 
ing <t; with the phonon-scattering conductivity a p h- The as- 
sumption of a homogeneous bec lattice with a low point- 
defect density gives a w a p h ■ Useful summaries of the prop- 
erties of a ph for such a system has been given by Urpin & 
Muslimov (1992) and by Pethick & Sarling (1995). Umk- 
lapp processes are by far the more important contribution 
to resistivity but become energetically disallowed at temper- 
atures T ~ Tu, where 

— «*w(if(i£)K 

and pi4 is the matter density in units of 10 14 g cm~ 3 . At 
temperatures Tu < T <ti Td, where Td is the Debye tem- 
perature, 



>ph 



5.5xl0 23 oL /6 fi^V /3 T Q 
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whereas in the limit TcTy, 

o ph = 2.1 xl0 28 p 8 f W3 T 9 



(19) 



(20) 



In the neutron-drip region p > p nd , the best value of the 
mass number A, on which the phonon spectrum depends, is 
not at all obvious. Pethick & Sarling suggest A — Aws but 
give no reasons why this choice should be suitable for super- 
fluid neutrons at T < T^ n . This factor introduces consider- 
able uncertainty in the form of the function a p h{T). It is also 
true that the modified phonon spectrum of an amorphous 
solid must produce substantial changes from equations (18) 
to (20). These have not been investigated here because, for 
the amorphous and Z-heterogenous solid, the impurity re- 
sistivity is large so that Oi < a ph except at temperatures 

T ~ 10 9 K. For comparison with equations (19) and (20), 
evaluation of equation (17) gives a conductivity 



<Ti = 2.0 x 10 2 



( Zpi 



\ Aws 



1/3 



(21) 



The values of Q given in Table 2 lead to high impurity re- 
sistivity which is only weakly-dependent on Aws- The fact 
that it is also temperature-independent removes a substan- 
tial degree of uncertainty from any model calculations made 
with it. 

Thus the (zero-field) electrical conductivity in most of 
the inner-crust volume is a w 10 s~ x at all times. This is 
particularly relevant to studies of the evolution of neutron- 
star magnetic fields which include both the Hall effect, with 
the possibility of a cascade to high wavenumber field com- 
ponents, and ohmic dissipation (see Hollerbach & Rudiger 
2002 and Geppert & Rheinhardt 2002, who also review ear- 
lier work on this problem). The Hall parameter considered 
by Hollerbach & Rudiger is 



-2-?- = 2.3 x 10 _3 Bi2 , - - - : 
n e ec \ ZpiA I 



f Aws\ 
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(22) 



12 P. B. Jones 



where n e is the electron density and B12 the magnetic flux 
density in units of 10 12 G. It exceeds unity only for B ~ 10 14 
G. The impurity resistivity found here is also relevant for any 
process in which the field distribution in the solid crust de- 
velops high-multipole components by mass-accretion or by 
flux expulsion from a superconducting neutron-star interior, 
or for models in which it is confined to the crust (see Page, 
Geppert & Zannias 2000, who also give an account of pre- 
vious work on distributions of this kind). Field components 
of wavelength 2h decay ohmically with an exponential time- 
constant 



which, for h — 10 5 cm, is only 4.5 x 10 5 yr, a time small 
compared with those usually considered in relation to decay 
of neutron-star internal or surface fields. Its temperature- 
independence would also simplify evolutionary calculations 
such as those of Page et al. For the reasons summarized 
here, we suggest that the structure of the solid phase of 
neutron-star matter is not an irrelevant detail. The standard 
assumption of a homogeneous bcc lattice is not adequate 
and the consequences of an amorphous and ^-heterogeneous 
structure should be included in very many different studies 
of neutron-star physics. 
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